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The Optical Variability of the Quasar 3C 279: The Signature of a 

Decelerating Jet? 

M. Bottcher and D. Principe^ 
ABSTRACT 

A recent optical monitoring campaign on the prominent quasar 3C 279 re- 
vealed at least one period of a remarkably clean exponential decay of monochro- 
matic (BVRI) fluxes with time, with a time constant of = 12.8 d, over about 
14 days. This is clearly too long to be associated with radiative cooling. Here 
we propose that this may be the signature of deceleration of the synchrotron 
emitting jet component. We develop a model analogous to the relativistic blast 
wave model for gamma-ray bursts, including radiative energy losses and radia- 
tion drag, to simulate the deceleration of a relativistically moving plasmoid in the 
moderately dense AGN environment. Synchrotron, SSC and external Compton 
emission are evaluated self-consistently. We show that the observed optical light 
curve decay can be successfully reproduced with this model. 

The decelerating plasmoid model predicts a delayed X-ray flare, about 2-3 
weeks after the onset of the quasi-exponential light curve decay in the optical. A 
robust prediction of this model, which can be tested with Fermi and simultaneous 
optical monitoring, is that the peak in the 7-ray light curve at ~ 100 MeV is 
expected to be delayed by a few days with respect to the onset of the optical 
decay, while the VHE 7-rays are expected to track the optical light curve closely 
with a delay of at most a few hours. 

Subject headings: galaxies: active — Quasars: individual (3C 279) — gamma- 
rays: theory — radiation mechanisms: non-thermal 
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Introduction 



The quasar 3C 279 {z = 0.536) is one of tlie best-observed fiat spectrum radio quasars, in 
part because of its prominent 7-ray flare shortly after the launch of the Compton Gamma-Ray 
Observatory (CGRO) in 1991. It has been persistently detected by the Energetic Gamma- 
ray Experiment Telescope (EGRET) on hoard CGRO each time it was observed, even in 
its very low quiescent states, e.g., in the wint er of 1992 - 1993, and is known to var y in 



7-ray flux by roughly two orders of magnitude (IMaraschi et al.lll994j : IWehrle et al.lll998l ). It 
has been monitored intensively at radio, optical, and more recently al so X-ray frequencies 



and has been the sub j ect of intensive mul tiwavelength campaigns (e.g.. IMaraschi et al.lll994 



Hartman et al.l 1 19961 : IWehrle et al.l Il998l ) . The most recent multiwavelength campaign on 
3C 279 included a Wh ole Earth Blazar Telescope (WEBT) campaign in the spring of 2006 
( iBottcher et al.l l2007bl ). During this campaign, the source was overall in a high optical 
state, with R ^ 14.0 - 14.5. However, the light curves showed an extraordinary feature: 
An unusually clean, quasi-exponential decay of the BVRI fluxes with a time scale of = 
12.8 d, extended over about 2 weeks. This paper aims at a theoretical interpretation of this 
extraordinary light curve feature. 

Flat-spectrum radio quasars (FSRQs) and BL Lac objects are active galactic nuclei 
(AGNs) commonly unifled in the class of blazars. They exhibit some of the most violent high- 
energy phenomena observed in AGNs to date. Their spectral energy distributions (SEDs) 
are characterized by non-thermal continuum spectra with a broad low-frequency component 
in the radio - UV or X-ray frequency range and a high-frequency component from X-rays to 
7-rays. In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in a 
relativistic jet. The high-frequency (X-ray - 7-ray) emission could either be produced via 
Compton upscattering of low frequency radiation by the same electron s responsible for the 
synchrotron emission (leptonic jet models; for a recent review see, e.g., iBottcherl l2007al ) . or 
due to hadronic processes initiated by relativistic p rotons co-accelerated with the electrons 



(hadronic models, for a recent discussion see, e.g., iMiicke fc Protherod I2OOII : iMiicke et al. 
20031). Several authors have modeled broadband SEDs of 3C 279 in variou s states (e.g., 
BednareklEggsl : Isikora et allboOll : Iflartman et aPboOll : iModerski et alll2003[ ). A consistent 
picture emerges that the X-ray - soft 7-ray portion of the SED might be dominated by 
synchrotron self-Compton (SSC) emission, while the EGRET emission might require an 
additional component, most likely external Compton emission. 

Standard leptonic models of blazar emission generally assume that a relativistic plas- 
moid containing ultrarelativistic nonthermal electrons moves with constant bulk Lorentz 
factor r along a jet, directed at a small angle with respect to our line of sight. How- 
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ever, for several blazars, in particular high-frequency peaked BL Lac objects detected at 
> 100 GeV 7-rays, such models sometimes require unexpectedly large bulk Lorentz factors 
(r > 50) and accordin g ly small viewing angles in order to explain their SEDs and variability 
teegelman et aliboosl : iGhisellini fc Tavecchiol boosi : IPinke et aliboosj ). Such large Lorentz 
factors and small viewing angles pose serious problems for AGN unification schemes, accord- 
ing to which FR I radio galaxies are believed to be the unbeamed equivalents of BL Lac 
objects. A possible solution to this dilern ma might lie in the deceleration of the emission 
region ( Georganopoulos fc Kazanas 2003a| jbl) from sub-pc scales, at which the optical - X- 
ray - 7-ray emission is produced, towards pc and kpc scales, which can be resolved with 
VLBA / VLB! techniques. At those scales, superluminal speeds of individual jet compo- 
nents of /5app ^ 10 are characteristically observed in most cases, providing an estimate of the 
Lorentz factor of jet components at those scales of F ~ 10. In fact, extreme deceleration of a 
radio-emitting plasmoid (component C3) in the jet of 3C 279 niay alr eady habe been directly 
observed in space VLBI monitoring observations (jPiner et al.ll2000l ). although the identifi- 
cation of this component over multiple observing epochs with different instruments/arrays 
is highly uncertain. 

In this paper, we propose a model analogous to the relativistic blast wave model 
which has successfully predicted and expl ained the smooth, self-si r ailar light curves of X- 
ray and optical afterglow s of 7-ray bursts (jPaczyhski fc Roadslll993l : iMeszaros fc Reeslll997l : 
Chiang &: Dermerlll999l ). We adapt this model for the specific situation in blazars. In par- 
ticular, we include self-consistently radiative losses and radiation drag from Comptonization 
of external radiation fields. A similar study, with emphasis on the details of the isotropiza- 
tion of particle distributions in the plasmoid an d on spectral feature s frorn various leptonic 
and hadronic processes, has been performed by iPohl fc Schlickeisen (l2000l ). who find good 
agreement of their results with characteristic SEDs of blazars. Here, we adopt a simpli- 
fied description of the particle dynamics and radiation processes, and focus on the expected 
monochromatic light curves dominated by the plasmoid deceleration. We review the ob- 
servational motivation from 3C 279 in ^ describe the model for the plasmoid dynamics 
in ^ and outline our treatment of radiation processes in §H As a test of our numerical 
simulations, we develop an analytical solution to the plasmoid dynamics and light curves in 
the self-similar deceleration phase in §3 In §2] we present results of our simulations and fits 
to the observed exponential flux decay of 3C 279 in January 2006. We summarize in ^ 

Throughout this paper, we refer to a as the energy spectral index, F,y [Jy] oc z/^". A 
cosmology with Qrn = 0.3, = 0.7, and Hq = 70 km s~^ Mpc~^ is used. In this cosmology, 
and using the redshift of z = 0.536, the luminosity distance of 3C 279 is = 3.08 Gpc. 
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2. Observational Motivation 

3C 279 was observed in a WEBT campaign at radio, near-IR, and optical frequencies, 
throughout the spring of 2006. Details of the observations, data analysis, and implications 
of the optical v ariabil ity patterns observed during that campaign have been published in 



Bottcher et al.l (l2007bl ). Fig. [T] shows the optical light curves of 3C 279 during spring 2006. 
The light curves exhibit an extraordinarily clean quasi-exponential decay with a character- 
istic time scale of ~ 12.8 days around JD 2453743 - JD 2453760. This light curve feature 
can not be interpreted as the signature of radiative cooling since the synchrotron cooling 
time scale for electrons emitting synchrotron radiation in the optical R band is 

r°y^^ ~3 X 10^5-'/'Di"'/'s (1) 

where Bq is the magnetic field in Gauss and Di is the Doppler factor in units of 10. This is 
of the order of at most a few hours for typical values of the magnetic field strength expected 
in quasars (5 ~ 1 G). Setting the synchrotron cooling time scale equal to the observed 
exponential decay time scale, would require a mag netic field of 5 ~ 7 x l^}-^ D^^^ G, which 
is about three orders of magnitude lower than usually inferred for quasar jets. We therefore 
favor a model in which the light curve decay is associated with the dynamics of the emission 
region rather than microscopic processes. We note that similar quasi-exp onential decays 



have also been observed in 3C 279 repeatedly in the 2007 observing season (ILarionov et al. 



20081). 



3. Model of a decelerating jet 

Our treatment of a decelera ting jet is borrowed froni the blast wave model of gamma- 



ray bursts. For details see, e.g., IChiang &: Dermerl (119991 ). We assume a plasmoid moving 



ballistically with initial mass Mq and bulk Lorentz factor Fq along the jet. Let M be the 
relativistic mass of the the plasmoid in the rest-frame of the plasmoid, then the momentum 
P of the plasmoid in the stationary AGN frame is P = /? F M c, where j3 is the normalized 
velocity v/c corresponding to the bulk Lorentz factor F. If the plasmoid radiates isotropically 
in its rest frame, the equation of motion of the plasmoid can be derived from momentum 
conservation, dP/dt = 0. However, in the case of a quasar, a substantial contribution to 
the (bolometrically dominant) 7-ray emission results from Compton upscattering of external 
radiation fields (EC = External Compton), and a significant transfer of plasmoid momentum 
to Compton-scattered external radiation ("Compton drag") has to be taken into account. 
We can therefore write 



dt \ dt J (3 

For large enough F, most of the EC radiation will be beamed into a narrow cone of solid 
angle Vt ~ 1/F^, and we can write the momentum transfer to EC radiation as 



If F2 . 

- / dL/dn cos edn^ — e'^c (3) 

c J Ati c 

where -E'ec is the internal energy loss due to EC radiation in the co-moving frame. For 
the purpose of an approximate, quantitative analysis to extract the salient spectral and light 
curve features of this model, we assume that all Compton scattering occurs in the Thomson 
regime so that 




^'ec = 3 car w'ext j N^ii) 7' d-^. (4) 

1 

The factor 1/F in Eq. (jlj) stems from the fact that E' constitutes a derivative with respect 
to time in the stationary AGN frame, and dt' = dt/T. We can also use Eq. ([3]) for a rough 
estimate of the magnitude of the radiation drag force, assuming that the observed 7-ray 
emission results from Compton scattering of an isotropic radiation field. Then, 



f) « -ii±£)|i:5?^ ^4,3X10- A, r- dyne (5) 

where /13 = uF^^'^'^ /{lO^^ Jy Hz) and Fi = F/10. This may be compared to an estimate 
for the Compton rocket effect due to radiation from the accretion disk. If we approxi- 
mate the accretion disk radiation as a point source (the most optimistic estimate) with 
luminosity Lo = L^q erg s~^, impinging from behind on an emission region of radius 
Rh = 10^^ -R16 cm, located at a distance r = 0.1 r_i pc from the accretion disk, the force on 
the plasmoid due to the Compton rocket effect can be estimated as 

(^) " " 2.3 X 10^^ Rl, rzl dyne. (6) 

Thus, for standard parameters with r > 10~^ pc and L^, < 10^^ erg s~^, the Compton rocket 
effect may safely be neglected in our calculations. 



The accumulation and radiative loss of relativistic mass M can be calculated as 
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M = A{r)pir)r{r) dr/dt + ^E',,^ (7) 

where dr/dt — /3c, A{r) = tt Rl{r) is the cross section of the jet, p(r) is the density of 
external material being swept up by the plasmoid, and 

oo 

E'rsd = ~ \caT u' J 7V,(7) 7' d^ (8) 
1 

Here, u' is the sum of the energy densities, u'q + Wg^t + ^sy) ^^d, again, we have assumed 
Compton scattering to be dominated in the Thomson regime. 

Derivatives with respect to time can be converted to derivatives with respect to distance 
r from the central engine, yielding 

dr M{r) dr At: M ^ ' 

and 

^ = ^«PWrM + ^(^ (10) 

We assume that a fraction of the swept-up electrons will be instantaneously accelerated 
into a power-law, which can be described by an injection function of the form 

Q(7) = Qo7~^-f^(7min,7,7m£«), (H) 

where H{xq,x,xi) — 1 ior xq < x < xi and otherwise. The low- and high-energy cutoffs 
of the electron injection function are given by 

where Se is the fraction of swept-up power that is transferred to rclativistic electrons, 
is the fraction of swept-up electrons which is accelerated to ultrarclativistic energies. The 
maximum Lorentz factor can be estimated by balancing the fastest conceivable acceleration 
time scale (the Larmor time scale) with the synchrotron loss time scale: 

7max~4.7x IO^^g'/' (13) 
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where Bq is the magnetic field in Gauss. The magnetic field may be parameterized in 
terms of a fraction of the swept-up energy transferred to magnetic-field energy density: 



B{r) = v/32 TT p(r) c2 T{r) (14) 

where an additional factor of 4 has been introduced to account for the compression of 
the pre-shock material by the strong shock. We can find the normalization Qq of the electron 
injection function through 



mm 



Eq. ([T]) indicates that the characteristic radiative cooling time scales of particles emit- 
ting synchrotron radiation at optical or higher frequencies are likely to be much shorter 
than the dynamical time scale of the system. Therefore, the balance of relativistic particle 
acceleration, cooling, and escape on a time scale 

rL = V ^ (16) 

yields a critical electron energy 7^ beyond which particles effectively radiate their energy 
away on a time scale shorter than the escape time scale: 

ic = -. r-B- (17) 

Aaru' rj Kb 

The resulting quasi-equilibrium electron energy distribution will be a broken powerlaw with 
parameters depending on whether 7c > 7min (slow-cooling regime) or 7c < 7min (fast-cooling 
regime). Re-writing the radiative cooling rate as 7 = — z/o7^ with z/q = (4/3) caxu' /{nieC^) 
we find for the slow-cooling regime 

7-9 for 7^in < 7 < 7c 

Nsch,t)^{ (18) 



Qo_ ^-(g+i) for 7c < 7 < 7n 



(9-1) 



In the fast-cooling regime, we have 



^j;;^) iLn 7 ^ for 7c < 7 < 7mi„ 
NU7,t)^{ (19) 



for7^i,<7<7 



max 



- 8 - 



The apparent large break in the spectral index around the synchrotron peak in the SED 
of 3C 279 indicates that this peak can not be attributed to a cooling break, but rather to a 



large value of 7miru and the system is in the fast-cooling regime (iBottcher. Reimer. fc Marscher 



20081 ). Eqs. ffT^ and ffT7|) naturally yield values that support the assumption of the fast- 



cooling regime for plausible parameter values: 



where Ce-i = Ce/O.l, while 



eOO^Fi (20) 



lc<SOTfR^,' (21) 



where we have only taken into account Compton cooling on the external radiation field for a 
characteristic value of u'^^^ as given in Eq. (l42l) . Consequently, 7min > 7c as long as (formally) 
r > 0.5 ^e/ce,-!, which will be the case for all plausible parameter values, even when the 
plasmoid becomes non-relativistic. 



4. Radiation 

The relativistic electrons described by the distribution functions (fT8l) and (fTOl) . will emit 
synchrotron, synchrotron-self- Compton (SSC) and EC radiation. For simplicity, we express 
the synchrotron emissivity using a 6 function approximation: 

.,,^.^^,.J^^4^(^"\JS\ (22) 



GirdleB {I + z) \eB J yVes 

where = 5/5„it with Ecnt = 4.414 x 10^^ G, and e' = ([1 + z]/D) hv^hs/ {^^0^). The 
corresponding photon number density is 

Q / / / \ 1/2 



Writing the electron distribution ( 1181) or ( 1191) as 



for 7i < 7 < 76 

iVe(7) = No { (24) 

for 7fe < 7 < 72 
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with (71,76,72) = (7mm, 7c, 7max) and {pi,p2) = {q,q + 1) for the sfow-coohng regime and 
(71,76,72) = (7c, 7mm, 7max) and (pi,P2) = (2, g + 1) for the fast-coohng regime, we can 
express the photon energy density as 



u 



- = ¥^^^^-^^1,^^ 3-pi 3-P2 



(25) 



The photon number density can be used in the IJoned (119681 ) formula to evaluate 
the SSC flux: 



obs J 



A-K dl{l + z) 



d^N,{^) / c/e'sy<y(e'sy)^?(e',e'sy.7) (26) 



with 



and 



^(e', e',y, 7) = ^ . f I — 1 I if -^7 < e < e' 



2 ^4p/ \ 



^(e', e'sy, 7) = ( 2 g In g + (1 + 2g) (1 - g) 



sy. 



/ fc sy 



472 
(4 6'sy7g)^ 



:i + 4eV7g) 2 



if 



e'sy < e' < 



4 6^7^ 
l + 4eV7 



where 



4eV7Ml- 



(27) 



(28) 
(29) 



We evaluate the external- Compton photon spectrum with a 6 function approximation 
for the external radiation field, 



n'ext(e, fi) 



u 



cxt 



■(5(e-e'ext)5(/i' + l) 



(30) 



where /i' = cos^'ext refers to the angle of incidence of the external photons with respect 
to the jet axis, and e'ext = Tcext- The Compton cross section is also approximated by a 5 
function, 



^ ar Sin' - n',) 5{e' - 7^ e'ext [1 - 
With these simplifications, the EC fiux can be calculated as 



(31) 
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5. Asymptotic Behaviour in the Deceleration Phase 

As in the well-known case of expanding blast waves in Gamma- Ray Bursts (GRBs), the 
plasmoid in a decelerating jet (with constant cross section A) is starting out in a coasting 
phase, in which the initial mass Mq greatly exceeds the swept-up relativistic mass in the co- 
moving frame. During this phase, the effect of the inertia of the swept-up mass is negligible, 
and the Lorentz factor remains roughly constant. This phase is followed by the deceleration 
phase. In the asymptotic limit of that phase, the initial mass of the plasmoid becomes 
negligible. 

Although we have properly included the effect of radiative cooling and radiation drag, in 
most cases the fraction of swept-up energy which is transferred to ultrarelativistic electrons 
and can therefore be radiated away efficiently, will be small. Therefore, we can approximate 
the equation of motion of the plasmoid by an adiabatic solution, as 



dr M 
and 



M(r)^Ap^ jT(r')dr' (34) 



ro 



where we have, for simplicity, assumed a constant external density p(r) = pext = const. 
This system has a self-similar solution of the form 



/ \ -1/2 

r(r) = To (^^ j . (35) 

Assuming for the purpose of an analytic estimate that we are looking right down the jet 
(^obs = 0), the observer's time as a function of Lorentz factor can be expressed as 



''obs 



r 



ro 



/3(r0 cos^obs) , [ dr' 

/3{r')c 2c J r2(r') 

ro 



4ro rgc' 



(36) 
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This yields a solution for the Lorentz factor as a function of observer's time: 

^0 ( '^o V^' 




The observed steep spectral index of the optical synchrotron emission from 3C 279 (a ~ 
1.7), indicates p ~ 4.4. This, in turn, signifies that the system is in the fast cooling regime 
since otherwise a cooling break would not produce a vFy peak at the synchrotron frequency 
corresponding to 7;,. Furthermore, electrons synchrotron radiating at optical frequencies 
are most likely beyond the break energy, i.e., 7 > 7min- For the prediction of synchrotron 
light curves we may therefore use Eq. [22] together with the lower branch of Eq. [191 For 
7max ^ 7min, the Coefficient Qo in Eq. [19] may be approximated as 

Qo(r) - ^-^^f;^^-"' (g - 1) l'^l{r) OC r^. (38) 

The cooling coefficient vq is expected to be dominated by synchrotron and / or external Comp- 
ton cooling. With the magnetic-field scaling from Eq. [HI both m^^^ and u'q carry a depen- 
dence OC r^. Thus, we find for the normalization of the ultrarelativistic particle population: 

iVe(7, W) OC 7-^^-^'^ OC tit'^)/^ 7-^^+'^. (39) 

In order to use Eq. [22]for a light curve estimate, we assume, again, for simplicity, ^^obs = 
and therefore D '^2T. Consequently, the characteristic electron energy 7 = \/e' /es oc V^^. 
This yields an expected light curve decay in the fast-cooling synchrotron regime of 

.F-(.ob., W) OC .S^)/^ (^-^) OC ^^/^ Kt^'^'- (40) 

In particular, for an injection index of g = 3.4, as inferred from the optical spectral index, a 
light curve of oc t~i^^ is expected. However, it should be pointed out that this can only be 
considered an upper limit to the steepness of the decay. Any non-zero observing angle will 
fiatten the decay of the light curve as it introduces a shallower decay of the Doppler factor D 
with decreasing F and therefore with time. In any case, by the time the plasmoid is in the self- 
similar deceleration phase, the monochromatic flux has already decreased by about an order 
of magnitude from its initial peak value, and is likely to be overwhelmed by other emission 
components in the jet. We do therefore not expect to observe the limiting deceleration 
case directly. We only developed this analytical case to demonstrate the agreement of our 
numerical simulations with the analytical expectation in the following section. 
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6. Numerical Results 

In order to highlight the salient features of our model, we assume a simple cylindrical 
jet geometry with a constant cross section of the jet, A = nRf, as well as a homogeneous 
external medium with density p(r) = pext- We achieved good fits to the observed optical 
light curves and overall SED shape of 3C 279 with the parameters listed in Table [H Our 
choice of the external density corresponds to a number density of next = 100 cm^^. This is 
at least about an order of magnitude lower than typical particle densities in the broad line 
regions of quasars. This is quite reasonable since we expect that the jet trajectory is already 
partially evacuated from previous ejection events. 

Fig. [2] compares snap-shot SEDs of our simulation to various observed broadband SEDs 
of 3C 279. It can be seen that the optical continuum spectra in various intensity states are 
well represented by the model, the X-ray spectral slope corresponds to the characteristically 
observed shape, and for much of the plasmoid evolution, the simulated X-ray flux is in the 
range of observed values. The 7-ray flux represented by our simulation corresponds to a 
medium to low state of 3C 279. Note that this work is not intended to attempt a model 
i nterpretation of t he very-high-energy 7-ray flux detected by MAGIC on Feb. 23, 2006 



(lAlbert et al.ll2008l ). in particular since we did not take into account any effects of intrinsic 



or intergalactic 77 absorption . For model implications of this high-energy detection see 



Bottcher. Reimer. fc Marscherl (120081 ). 



The upper panel in figure [3] illustrates the evolution of the plasmoid bulk Lorentz factor 
r as a function of distance from the central engine, obtained by numerically solving the 
coupled system of Eqs. and ffTU]) . 

As in the well-known case of the blast-wave model for gamma-ray bursts, the plasmoid 
evolves through a coasting phase with approximately constant Lorentz factor Fq. Around 
the characteristic deceleration radius r^, at with the swept-up relativistic mass equals the 
mass of the initial ejecta Mq, the evolution makes a gradual transition into the asymptotic 
self-similar deceleration phase with F(r) oc r~^/^, treated analytically in §3 For the input 
parameters listed in Table [1], this occurs at 

rd = ~ 1.1 X 10^^ cm. (41) 

LqIT Kf^ Pext 

Note, however, that substantial deceleration happens long before is reached. The 
lower panel of Fig. [3] shows the resulting R-band light curve over the entire evolution of the 
plasmoid, down to a mildly relativistic speed. Again, we see the gradual turnover to the 
self-similar deceleration phase with the expected light curve decay as Fgy oc t~^^'^ (see 
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During the transition from the coasting to the self-similar deceleration phase, the light curve 
is reasonably well represented by an exponential decay. 

The detailed fits to the V, R, and I band light curves of 3C 279 during the quasi- 
exponential decay around January 15, 2006, are shown in Fig. |H The B and U band light 
curves were rather poorly sampled during this time period so that the (equally good) fits to 
those light curves does not provide substantial additional information. 

Using characteristic parameters for the accretion disk luminosity, = 10^^ L45 ergs s~^ 
and the Thomson depth tt^blr = 10^^ r_i and radius -Rblr = 0.1i?_i pc of the BLR of 
3C 279, we estimate an external radiation energy density of 



Wext ~ .^""p™ ~ 3 X 10-3 L45 r_i RZI ergs cm'^ (42) 

4 -Kblj^ C 



and ecxt ~ 10 ^ (see, e.g., iBottcher. Reimer. &: Marscherl |2008| ). which we use for our cal- 



culation of the EC 7-ray emission component and the associated radiation drag term. The 
predicted light curves of the combined SSC + EC emissions at X-rays and 7-rays are com- 
pared to the R-band light curves in Fig. [51 

The most remarkable feature of the predicted light curves is a delayed X-ray outburst 
about 2-3 weeks after the onset of the optical decline. This is a consequence of the shift 
of the SSC peak towards lower frequencies as the blob decelerates, which follows from the 
dependence t'pk.ssc oc j^^^ oc F®. This SSC peak frequency decrease is much more rapid 
than the decrease of the SSC peak fiux. This spectral evolutio n of the SSC component is illus- 



trate d in Fig. [2J In fact, RXTE monitoring of the source (see lBottcher. Reimer. &: Marscher 



20081 ) did detect a transition of 3C 279 from a quiescent state throughout the first half 
of January 2006, to a very active, high X ray fiux state with a substantial X-ray outburst 
around Feb. 5, 2006. It is tempting to postulate that this may have been the X-ray signature 
of the plasmoid deceleration observed during the optical decay around January 15. 

A prediction of this scenario is that the ~ 100 MeV 7-ray light curve peaks a few days 
after the optical before it begins its quasi-exponential decay. This is a consequence of the 
fact that this energy range consists of comparable contributions from the SSC and the EC 
radiation components. As the SSC component shifts towards lower energies during the early 
deceleration, the EC component just moves into the 100 MeV range, before it also decays 
away to lower peak energies and lower fiux levels. The very-high-energy (> 100 GeV) light 
curve is expected to follow closely the optical one, with only a small delay, which critically 
depends on the initial bulk Lorentz factor and the mean photon energy of the external photon 
field. 
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We caution that any predictions concerning absolute flux levels at X-rays and 7-rays 
are very model-parameter dependent and vary substantially for different choices of initial 
mass, initial bulk Lorentz factor, initial pladmoid radius, external matter density, external 
radiation field, etc. Furthermore, we have not in cluded any effects of intrinsic (in partic- 
ular, in the radiation fiel d of the BLR, see, e.g., iDonea fc Protherod l2003l : iReimen 120071 : 



Sitarek fc Bednarekl l2008l ). nor intergalactic 77 absorption. However, while this may sig- 



nificantly affect the overall 7-ray flux level, the light curve features discussed above are 
predominantly a consequence of the plasmoid dynamics, which are dictated by the observed 
optical light curves, and are therefore robust predictions of the decelerating-jet model. 



7. Summary and Conclusions 

Motivated by an extraordinarily clean quasi-exponential decay of the V, R, and I band 
light curves of 3C 279 over a period of about 2 weeks during a recent WEBT campaign, we 
proposed a model of a decelerating plasmoid in the jet of this quasar. We take into account 
self-consistently the inertia of swept-up mass as the plasmoid propagates through the gas of 
the AGN environment, as well as radiation drag and radiative cooling. We have demonstrated 
that, similar to the relativistic blast wave model for GRBs, the plasmoid makes a transition 
from a coasting phase with approximately constant Lorentz factor, to a self-similar phase. In 
the case of a homogeneous external medium and a cylindrical jet, the self-similar deceleration 
phase is described by r(r) oc r~^/^ oc tj^l^. The resulting optical synchrotron light curves 
are well approximated by quasi-exponential decays during the transition from the coasting 
to the self-similar decay phase. In the asymptotic limit of the self-similar decay phase, the 
synchrotron light curves at frequencies corresponding to electron energies above the cooling 
break of the electron spectrum follow a behaviour of z^-F;y(t'obs i^obs) oc t'obl'^''^ '^obs^^'^^^^; where 
q is the injection spectral index of ultrarelativistic electrons in the plasmoid. 

We note that the choice of a conical jet instead of a cylindrical one would recover the 



isotropic blast wave model for gamma-ray bursts (e.g., IChiang &: Dermerl Il999l ): It would 
lead to a scaling of the plasmoid surface as A oc r^. In the self-similar deceleration phase, 
this recovers the well-known adiabatic blast wave solution with r(r) oc r~^/^ oc t^^l'^ . This 
case describes a much faster deceleration of the plasmoid and therefore substantially steeper 
light curves than the cylindrical jet geometry assumed here. 

We have demonstrated that this model can adequately reproduce the observed optical 
light curves of 3C 279 during the ~ 2 week long quasi-exponential decay phase in January 
2006. This model predicts a delayed (SSC-dominated) X-ray outburst about 2-3 weeks 
after the onset of the optical decay. We speculate that the X-ray flare around Feb. 5, 
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2006, detected by RXTE monitoring, may have been the X-ray signature of the plasmoid 
deceleration seen earher in the optical bands. 

A robust prediction of the decelerating plasmoid model which can be tested with Fermi 
and simultaneous optical monitoring is that the peak in the 7-ray light curve at ~ 100 MeV 
is expected to be delayed by a few days with respect to the onset of a quasi-exponential 
light curve decay in the optical, while the VHE 7-rays are expected to track the optical light 
curve closely. 

The quasar 3C 279 is one of the most active blazars known. The ejection of 7-ray 
emitting plasmoids from the nucleus of 3C 279 might be a frequent event. We do therefore 
expect to observe the quasi-exponential decay phase following such ejection events only in 
a temporary quiescent phase, in which it is not overwhelmed by subsequent ejection events. 
This might be the reason why these quasi-exponential decays in the light curves of this blazar 
are not more frequently observed. In addition to one or two occurrences during the 2006 
WEBT campaign, several more quasi-exponenti ally decaying light cu rve segments can also 



be identified in the 2007 WEBT campaign data (ILarionov et al.ll2008l ). indicating that these 
are rare, but not unique events. 
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Fig. 1. — Light curves of 3C 279 in various optical bands during the spring of 2006. The 
dotted red hne indicates our model fit, as discussed in the following sections, to the R-band 
light curve during the quasi-exponential decay around January 15, 2006 
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Fig. 2. — Simulated snap-shot SEDs from our decelerating plasmoid model, compared to 
various observed SEDs of 3C 279 in various observing epochs and activity states. The time 
sequence goes from red - dotted — *• dashed — >• dot-dashed long-dashed — solid orange 
- dotted — > dashed — > dot-dashed long-dashed. This simulation provides a good fit to the 
optical light curve decay of 3C 279 around January 15, 2006 (see Fig. HI). 
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Fig. 3. — Top: Evolution of Lorentz factor F as a function of distance r from the central 
engine for the simulation illustrated in Fig. [2l Bottom: R-band light curve from this 
simulation. 
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Fig. 4. — Fit of our decelerating-jet model to the V, R, and I band fluxes of 3C 279. 



Table 1. Parameters of our plasmoid evolution simulation providing a fit to the V, R, and 
1 band light curves of 3C 279 around January 15, 2006. 



Parameter 


Symbol 


Value 


Initial Lorentz factor 


To 


50 


External matter density 


rioxt = Poxt / rrip 


100 cm-3 


Plasmoid radius 


Rb 


3 X 10^'' cm 


Initial mass 


Mo 


2.6 X 10^1 g 


Electron injection index 


q 


3.4 


B-field equipartition parameter 




10-3 


Electron acceleration efficiency 


ee 


0.1 


Relativistic electron fraction 


Ce 


0.5 


Electron escape time scale parameter 


V 


10 


External radiation energy density 


ttext 


3 X 10-3 erg cm-3 


Observing angle 


^obs 


5° 
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Fig. 5. — Predicted X-ray and 7-ray light curves for parameters from our fit to the optical 
light curves of 3C 279. 



